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Density waves of granular material (sand) flowing through 
a vertical pipe have been investigated. Clear density waves 
emerge when the cock attached to bottom end of the pipe is 
closed. The FFT power spectra were found to show a sta- 
ble power-law form P(f) ~ f~ a - The value of the exponent 
was evaluated as a = 4/3. We also introduce a simple one- 
dimensional model which reproduces a = 4/3 from both sim- 
ulation and theoretical analysis. 



P(f) ~ f~ a ■ The result is a = 4/3. We also propose a 
one-dimensional model supplemented by the white noise 
which reproduces a = 4/3 near the neutral curve of the 
linear stability analysis of uniform states. In other words, 
we will clarify the origin of the power-law in density waves 
of granular pipe flows. 
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Recently much attention has been paid to the dynam- 
ics and statistics of granular materials because of their 
ubiquity in nature and the application to technology. Un- 
like usual solids, liquids or gases, granular materials are 
known to show complex dynamical behaviors Jt| , such as 
convection 0, size segregation ||, bubbling [pj7 standing 
waves and localized excitations under vertical vibration 
HP] and a fluidized bed due to air injected inside a box 
containing granules f^M . 

Pattern formation of grains flowing through a vertical 
pipe which can be regarded as a one dimensional realiza- 
tion of a fluidized bed is also a typical example of unusual 
features of granular motion 1^13]] . Emergence of density 
waves (e.g. slugging) has been investigated by molecular 
dynamics (MP) and lattice-gas automata (LGA) simu- 
lations [^],[lO| and by the experiments using sand in air 
pd| , ^2[ and metallic spheres in liquids Power- law 

form of the power spectrum P(f) ~ f a , where / is fre- 
quency, of density fluctuations was also found in both 
numerical simulations || and experiments |Tl|JT^ ]. Al- 
though their interpretations on the origin of the emer- 
gence of density waves are different, estimated values of 
the exponent a is close to each other (1.3 < a < 1.5). 

Although the previous experiment |ll| , |l2| reported a = 
1.5, some of their experimental procedures seem a little 
ambiguous: Since they merely plugged up the bottom 
hole of a pipe by half in order to induce density waves, the 
rate of air flow out of the bottom end of the pipe was not 
well controlled. Besides, the power spectra they obtained 
were still noisy. In this Letter we will first present better- 
controlled air flow out of the pipe and more accurate 
experimental results than the previous ones by increasing 
the number of trials. One of our results is the precise 
estimation of the scaling exponent of the power spectrum 



FIG. 1. Schematic illustration of the experimental setup. 



Our experimental setup is shown in Fig. 1. We used 
a glass pipe of 1500 mm length and 3 mm inner diam- 
eter, i.e., the same one used in Ref. p] , ^2| . A flask is 
connected to the bottom end of the vertical pipe and 
a flow-meter is attached to an outlet stuck out of the 
flask. We used rough sand, i.e., our granular particles 
are polydisperse and their average diameter is about 0.3 
mm. Thus the inner diameter of the pipe is about ten 
times larger than the average diameter of sand. We pored 
rough sand into a hopper connected to the top end of 
the pipe, which flows through the pipe due to gravity. 
Sand finally falls into the flask while air can exit out of 
it through the flow-meter, which can control the rate of 
air discharge. Meanwhile we measured density fluctua- 
tions as the transmission light intensity across the pipe 
at a fixed location midway up the pipe by using laser 
light and detecting system (KEYENCE, LX2-02). The 
laser light has a rectangular cross-section of 10 mm wide 
and 1 mm high and is emitted with a pulse frequency of 
4096 Hz. The detecting system has linear response be- 
tween the output voltage and the density of particles in 
the range of output voltage between 1 and 5 volts. 

Let us now consider what happens by closing the cock 
of flow-meter. Granular particles falling in a pipe intcr- 
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act with the medium (air) due to the viscosity. When 
the flow-meter is removed from the apparatus, i.e., the 
bottom end of the pipe is fully open, granules can fall 
rather freely as if the existence of air could be neglected. 
(Hereafter we will refer this situation to fully open.) In 
fully open there are no visible density waves, as reported 
in Refs. Jll],[l2|]. The reason is that both granules and 
air flow together through the pipe. As the cock is grad- 
ually closed, however, the pressure in the flask rises and 
the effect of the viscous force becomes more important. 
In particular, when the cock is fully closed (we refer to 
fully closed), air in the pipe must go upward due to ap- 
proximate conservation of the total volume (sand plus 
air) in the flask while sand falls downward. We consider 
that the increment of interaction between granules and 
air induces density waves. 



the medium air was controlled to flow downward by grad- 
ually opening the cock of the flow-meter, the position of 
the emergence of clear density waves gradually shifted 
downward. The details of the results will be presented 
elsewhere.) 



FIG. 3. Log-log plot of power spectra P(f) of time series 
signals in fully closed. The straight line with the slope of —4/3 
in the figure is for an eye-guide. 



FIG. 2. Time series signal of the granular flows for fully 
open (top) and fully closed (bottom). Upper signal is shifted 
by 0.4 V to avoid the data overlap 

Figure 2 shows typical time series signals of granular 
flows measured at x — 100 cm, where x is defined as the 
distance of measuring point along the pipe from the top. 
In contrast with the white-noise-like signal (top) in fully 
open, one can see that the density wave for fully closed 
(bottom) has intermittent structure. In this figure higher 
(lower) voltage corresponds to smaller (larger) granular 
density since we measured the light intensity transmit- 
ted across the pipe. In particular, the output voltage is 
about 2.15 volts when a pipe is fully packed with sand, 
and about 2.50 volts when it is completely empty. This 
assures us that our measurements were well in the linear 
range between the voltage and the density. In the case 
of real flow patterns one can observe that two clusters 
sometimes collide with each other and merge into one, 
and sometimes one cluster splits into more than one clus- 
ter, in a behavior reminiscent of a chain of traffic jams on 
a crowded highway In the following discussion 

of this Letter we restrict ourselves to fully closed. (When 



Let us pay our attention to the power spectrum P(f) 
of density fluctuations of granular flows. Figure 3 shows 
P(f) of recorded signals in fully closed. Each spectrum 
was obtained by averaging over 640 independent data 
with length 8192 discrete points each (two seconds in real 
time) and was appropriately shifted to avoid the overlap. 
The self-organization of the power-law form P(f) ~ f~ a 
is observed as the measuring position x increases: The 
spectrum at a; = 10 cm is similar to that in fully open. 
Density waves have already emerged at x — 40 cm where 
the scaling regime is covered by the whole frequency 
range except for the fast decay in higher /. Then the 
system falls into the steady state from about x = 50 cm 
downward. The scaling range for the steady state, i.e., 
for x > 50 cm is from 10 Hz to 200 Hz, as seen in Fig. 3. 

In the frequency range from 15 Hz to 150 Hz, values of 
the slope are fitted by the least mean square method as 
a = 1.33 ± 0.04, 1.31 ± 0.05 and 1.34 ± 0.05 for x = 80, 
90 and 100 cm, respectively. (The data at x =90 cm and 
100 cm are not shown in Fig. 3.) These values are very 
close to 4/3 suggested by LGA simulation 

Let us allocate the rest of this Letter to introduce a 
one-dimensional model which reproduces a — 4/3 in fully 
closed system. 

The important mechanisms for particle dynamics are 
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the drag between air and particles, and the relaxation 
process to an optimal velocity which may be the sed- 
imentation rate of particles. When the N particles are 
confined in a quasi-one dimensional container, the motion 
of particles may be described by the following nonlinear 
equation: 



r n + C[r„ ~ W({r n })} = TF({r n }) + /„(*), 



(1) 



where r n and £ are the relative distance between the 
n-th and (n + l)-th particles, and the drag coefficient, 
respectively. The collisional force F({r n }) — p'(r n+ i) + 
<p'{r n -i)— ^• i p'{ r n) comes from a soft core repulsive poten- 
tial <p(r n ). The parameter T represents the strength of re- 
pulsion. The optimal velocity ^({rn}) = U( Tn+1 9 +,n ) — 



r n +r n 



— ) is the linear combination of sedimentation 



rate U(x) which is the nonlinear function of the local 
volume fraction [JL7| in general. The most crucial simpli- 
fication of (|l|) without any justification is to assume that 
/„ is a Gaussian white noise with zero mean. It should be 
noticed that the drag £ is irrelevant in fully open, because 
air in the pipe flows away together with particles. Thus, 
the motion of particles are almost elastic and clear den- 
sity waves cannot be observed in fully open. It should 
be noted that (Q) is written for the relative motion of 
particles. Discarding the noise term from (^), we obtain 
an equation of motion for particle at the position x n as 
x n + C[x n - U C n+1 l an - 1 )] = T[<p'(r n ) - <p'(r n -i)] where 

Linearizing ([!]) around the uniform solution r n = 
where a = f n = TV" 1 ~^2 n r n , we obtain 

r k + C[r'k - iU'shxkr k ] = 2Tcp"(cosk- l)f* + f k {t), (2) 

where the argument of U' and ip" is a. fu and fk(t) are 
respectively the Fourier transform of 5r n = r n — a and 
f n (t). Equation (j|) has the solution r k (t) oc exp[<j±t] , 
where 




-2Tip"(l-cosk)+i(U'smk. (3) 



Re[<7+] represents the relevant eigenvalue of the linear 
problem, which becomes positive for U'(a) 2 cos 2 (fc/2) > 
Tip" (a). Thus the most unstable wave number is /c — > 
and the neutral curve is given by T c = U' 2 /ip"(a). At 
T = T c (l — fi) the expansion of cr + around k — is given 

by a + {k) ~ l [c«k + ■■■] + ^fc 2 - ifc 4 + • • 

6 C _ 4C 

where cq = U'{a). Thus, for fj, > the uniform state is 
unstable due to the negative diffusion. 

Adopting U(r) = tanh(r— 2)+tanh(2), tp(r) = sech 2 (r) 
, £ = 2, N = 256, T c = 3.95798 • • •, and a = 2 at t = 0, we 
simulated (|l]) by the classical Runge-Kutta method until 
t = 2 11 with time interval At = 1/2 4 under the periodic 
boundary condition. We used the uniform random num- 
ber distributed between —X and X with X = 9/1024 
fist for f n (t). Figure 4 displays the power spectrum 



P(f) = |c(/)| 2 obtained from our simulation of (|T|) at 
/i = 1/64, where c(/) is the Fourier transform of the dis- 
cretely sampled data of the density c(t) = — ^ — ^— ■ 

with the interval 1. This clearly supports P(f ) ~ f~ 4 ^ 3 
as in our experiment. From the examinations of several 
values of ji, we have confirmed that the qualitative re- 
sults are insensitive to the sign of fi when \fi\ -C 1. This 
result is reasonable because near the neutral curve the 
time scale of relaxation or growth of fluctuations is much 
longer than the time scale induced by the noise f n (t). 
Our result suggests that the linear relaxation theory of 
fluctuations can be used to explain P(f) ~ f~ 4 ^ 3 - 



FIG. 4. Log-log plot of power spectrum P(f) obtained from 
the numerical integration of (1), where the unit of / is l/(27r). 
The guide line represents /~ 4//3 . See the text for details. 



Let us briefly sketch how to derive 4/3-law from the 
behavior of structure factor 



Sk{t) =J2< exp[ik{Sr n {t) - Sr m {0))] > 



(4) 



in weakly stable states, i.e., fi < and \/j,\ <C 1, 
where 8r n — r n — a is the fluctuation of relative dis- 
tance. The structure factor can be rewritten as Sk(t) = 



exp 



2 



where <j> nm (t) =< (Sr n (t) - 

Sr m (0)) 2 >. For fi < 0, 5fe(t) can be calculated as in 
the case of polymer dynamics |l9") . With the aid of the 
expansion of a + , (|l|) is reduced to 

d T r(z, t) - d 3 z r(z, t) = e[d 2 - 8*]r(z, r) + £(z, r) (5) 



where r = e 3 /?t, z 
with e = 



^e(x + Cot), £(z,t) = e 3 /3/„(t) 
3c 




[i and (3 



The solution of 



d) is given by f fc (r) ~ JJ" ds exp[Afc(r - s)]| fc (s), where 
Afe = ik 3 — efc 2 (l + fc 2 ). Thus, we obtain the correlation 
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< f fc (r)f_fc(0) >= 



D 



2elk 2 (l + k 2 ) 



exp[A 



(6) 



where I is the system size, and we used < Cfc( T )?p( r ') >— 

A long but straightforward calculation parallel to Ref. 
@ yields 



5fc(r) ~ y da; exp[— Dok 2 T 



DP 
~2e 



Dk 2 



.1/3 



h{u% 
(7) 



where the argument of Sk is replaced by the scaled time, 
Dq is the diffusion constant for the gravitational center 
in (^]), u — xr -1 / 3 , and w = \z — z'\. Since h(u) converges 
to h(Q) — 7r/r(l/3) as time goes on, we obtain 



S*(r) 



Dk' 



■exp[- 



Dk 2 



er(l/3) 



,1/31 



(8) 



in intermediate time range. In the limit of small r, 
ADk 2 

SkM oc 1 ^^^fc 2 T 1/3 H . Thus its Fourier trans- 

v 7 eT(l/3) 

form, giving the power spectrum Pk{f) obeys 



Pkif) ~ r 



4/3 (as /->«>), 



(9) 



where use was made of dte ,27r ^*|t| 1 / 3 oc / 4 / 3 . The 
value 4/3 is identical to the one obtained by our experi- 
ments and numerical simulations || . Thus our model (|) 
reproduces a = 4/3. It should be noted that the appear- 
ance of this power-law form in the original model (|l|) is 
only for / < £ since we eliminate the fast decaying mode 
cr_ in our analysis. This tendency is also observed as the 
higher- frequency cutoff in our experiment (see Fig. 3). 

The 4/3-law is determined by short time behavior of 
the dynamics of density waves induced by the noise, 
which is essentially determined by the linear dispersion 
relation Afc ~ ik 3 . The details of our theoretical analy- 
sis including the effects of nonlinearity will be discussed 
elsewhere pcfl . 

In this Letter we have confirmed a = 4/3 as the power- 
law exponent in the frequency spectrum of density corre- 
lation function from the experiment, the simulation and 
the theory. We have also clarified the mechanism to yield 
/~ 4 / 3 spectrum which is related to the critical slowing 
down of the density fluctuations. It should be noticed 
that the continuous increase of a in LGA jll| from a = 
to 2 with the particle density is consistent with the 4/3- 
law and our picture, because the spectrum determined 
by the noise in linearly stable uniform state far from the 
neutral curve should be white (a = 0) and the effective 
exponent of the power-law becomes large when the expo- 
nential decay (i.e. a = 2) in the off-critical region exists. 
There is, however, a discrepancy between our results and 
the one on experiments in liquids |b| . The reason of this 
difference should be clarified in the future. 
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